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A method is presented for deriving the nonrelativistic quantum hamiltonian of a free massive 
fermion from the relativistic lagrangian of the Lorentz-violating standard-model extension. It per- 
mits the extraction of terms at arbitrary order in a Foldy-Wouthuysen expansion in inverse powers 
of the mass. The quantum particle hamiltonian is obtained and its nonrelativistic limit is given 
explicitly to third order. 



I. INTRODUCTION 

Os . 
ON' 

■ Establishing the physical relevance of a lagrangian in relativistic quantum field theory often requires a determination 
of its nonrelativistic content. The Foldy-Wouthuysen (FW) transformation (!]] provides a systematic approach to 
O .■ understanding the low-energy effects of certain theories. Given the relativistic quantum hamiltonian for a theory 
of massive four-component fermions, the nonrelativistic quantum hamiltonian for the corresponding two-component 

m : 

particle can be derived in an expansion in inverse powers of the fermion mass. 

In this work, we use generalized FW methods to investigate the quantum particle hamiltonian that describes the 
physics of a free massive two-component fermion emerging from the relativistic lagrangian of the Lorentz-violating 
standard-model extension [0 . This standard- model extension is based on the idea of spontaneous Lorentz breaking 
in an underlying theory f| and has been used for various investigations placing constraints on possible violations of 
| Lorentz symmetry |^],|]-[l5[ , several of which depend crucially on the nonrelativistic physics of free massive fermions. 
^J - ■ In these investigations, specific terms in the nonrelativistic hamiltonian have been derived as needed, but a full 
treatment has been lacking. Here, we provide a systematic approach that permits extraction of the relevant terms 
in the nonrelativistic hamiltonian at arbitrary order in the FW approximation. We obtain the quantum particle 
hamiltonian and provide explicitly the form of the nonrelativistic hamiltonian to third order. Our results are directly 
. relevant to recent analyses of muon and clock-comparison experiments Jl^ , p^[ and are expected to have substantial 
""^ impact on further studies of the physical implications of the standard-model extension. 

The general form of the relativistic lagrangian for a free spin-^ Dirac fermion ip of mass m in the standard- model 
I ' extension is 

OS. 

C = !n&(7„ + + d^'f + e v + if vlb + \gx^<J Xtl ) d v ip - ^(m + + + ±H^a^)ip . (1) 



X 



C3 



This is a generalization of the usual relativistic lagrangian for a free massive Dirac fermion. The Dirac matrices 
{1, 75, 7 M , 757^, a^} have conventional properties, and the Minkowski metric r]^ has signature —2. The parameters 
a^, <7,„, d^v, e M , f^, gx^v, and control the extent of Lorentz violation in the theory. In a given observer inertial 
frame, they can be regarded as fixed real Lorentz vectors or tensors. Note that H^ v can be taken as antisymmetric, 
<7t„ and d^ v as traceless, and g\^ v as antisymmetric in the first two indices. Since Lorentz symmetry is known to 
be valid to high precision, any nonzero parameters in nature would need to be minuscule. We therefore restrict our 
attention in this work to terms linear in these parameters. 

In the next section, the relativistic particle-antiparticle hamiltonian H corresponding to the lagrangian ([!]) is 



obtained. Some basic information about our procedure for extracting its FW form is discussed in section III, together 
with our definition of the relevant FW sequence. Features of this sequence are derived in section |vL and the quantum 
particle hamiltonian and its nonrelativistic limit to third order are explicitly presented in section |V|. 



II. RELATIVISTIC QUANTUM HAMILTONIAN 



The first step in deriving low-energy effects of the Lorentz-violating terms is to obtain the relativistic hamiltonian 
H associated with the lagrangian ([!]). However, methods for direct construction of H are inadequate because Eq. (|l|) 
contains couplings involving time derivatives. For example, applying the Euler-Lagrange equations to £ and solving 
for H from the equation id^rp = Hip results in a non-hermitian hamiltonian and a corresponding nonunitary time 
evolution. 
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One method of bypassing this technical difficulty is to perform a field redefinition ip = A\ in the lagrangian [EJ, 
with A chosen such that the dependence of the lagrangian on do\ is just that of the usual Dirac lagrangian. Then, the 
wave function associated with \ evolves conventionally in time. The field redefinition leaves unchanged the physics, 
while it causes the time-derivative couplings to be replaced by extra terms in the lagrangian. 

To implement this idea, we write the lagrangian ([!]) in the forms 

C = \ii>Y v d v ip - ~$Mij) 

= yXJ d°x+k i x(Ar j A)& x-x(AMA) X , (2) 

where r„ and M are defined according to the correspondence with Eq. (Q), and ip = \A with A := j°A^°. In 
the second expression the Lorentz indices are separated into timelike and spacelike cartesian components, /i = and 
j = 1,2,3, with summation on repeated indices understood. 
The choice 

A = i-i 7 °(r - 7o ) , a = i- ±(r - 7o ) 7 ° (3) 

implements the equality (||) to linear order in the parameters for Lorentz violation. Derivation of the relativistic 
hamiltonian H can then proceed through the Euler-Lagrange equations, which take the form of a modified Dirac 
equation: 

{iAT^Ad^ - AMA)x = . (4) 

We find 

H = -fAT 3 Ap> + j°AMA , (5) 

where the three-momentum of the particle is denoted pj, and H obeys the equation idox — Hx- 
Explicitly, the relativistic hamiltonian can be written 



H = m( 7 ° + V a + O + So) , (6) 



where 



mV - - 
mSo 



= -PjH 7 

= [-b + (d j + djo^Ts + [ a j ~ (cjk - coorijk)p k ]l Q l } + ifjp'lhl + i[H j + (g 3 ok + 9jko)p k ]l J , 
= [ao - (coj + c j0 )p> - me ] + [-bj + (d jk - d 00 r] jk )p k - \me klm -q jm 9kio\lbl l 3 

-[mcoo + ej-p J ] 7 ° - [±e klm r)j m H k i + md j0 - e lmn r)j n {\gi mk - Vk m gm)p k ]l5l :i ■ (7) 

In these expressions, the totally antisymmetric rotation tensor e jW satisfies £123 = +1 and e^ kl = —Sjki- The particular 
decomposition of H into the four terms in (^|) is chosen for later convenience. 

As an aside, we remark that the relativistic hamiltonian is also readily found if the theory (Q) is extended to include 
a minimal coupling to a U(l) gauge field A^, It suffices to replace the partial derivative id^ in Eq. ([!]) with the 
covariant derivative iD^ := id^ — qA^, where q is the particle charge. The relativistic hamiltonian then has the same 
form as in Eqs. (|^) and (R) above, except that all occurrences of pj must be replaced with ttj :— pj — qAj and the term 
qAo must be added to Eq. (|^). The resulting hamiltonian is relevant, for example, for studies of Lorentz- violating 
effects in quantum electrodynamics. 



III. DEFINITION OF THE FW SEQUENCE 

In the strict nonrelativistic limit, the lower two components of the relativistic wave function x are negligible, so the 
upper two components of x suffice to determine the nonrelativistic particle behavior. However, more generally the 
Dirac equation couples the upper and lower components of x- The object of the Foldy-Wouthuysen procedure is to 
find a (momentum-dependent) unitary transformation 

H^H := e iS He~ iS = exp[a,d(iS)]H , (8) 
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where a.d(X)Y := [X, Y] , such that H is 2 x 2 block diagonal. This therefore decouples the upper and lower components 
of the FW-transformed wave function <f> := e lS \- Requiring hermiticity of S ensures that e lS is unitary. It follows that 
H is hermitian and that both hamiltonians H and H describe the same physics. The FW transformation amounts 
to a unitary rotation in the Hilbert space of the free particle states that preserves the dominance of the upper two 
components of the wave function. The quantum particle hamiltonian h re \ and the nonrelativistic limit h we seek are 
given by the leading 2x2 block of H . 

Solving directly for H would be of interest but is challenging in the general case. Instead, we present a method 
that allows approximation of H to arbitrary accuracy in an expansion in powers of \p\/m. The basic idea is to apply 
a succession of transformations of the type (||) , chosen so that each iteration of the transformed hamiltonian has a 
smaller block off-diagonal part than the previous one. The exact FW transformation is the limit of this sequence. 
Although more direct approaches can yield a low-order approximation to h without the use of our method, the results 
derived here permit straightforward calculation of h le \ and of h to any desired order. 

For definiteness in what follows, we work within the Dirac-Pauli representation of the Dirac matrices, for which 

7° = ^ q ^ ^ and 7 J ' = f ^ j °qJi wnere cr ' J are the usual Pauli matrices. We define a matrix to be even if it is 

block diagonal and odd if it is block off-diagonal. Any 4x4 matrix X can be uniquely written as the sum of an even 
part and an odd part, X = Even(A) + Odd(A), where Odd(A) = ±7°[7°, X] and Even(A) = ^7 {7°, X}. 

We seek a sequence of FW transformations such that the odd part of the hamiltonian progressively decreases in 
some suitable matrix norm, such as \\A\\ := max a .b{|^4 a b|} for a.b — 1,2,3,4. In the remainder of this section, an 
appropriate sequence {H n } of hamiltonians is introduced. For each n, we also introduce a parameter t n that turns 
out to provide a measure of the size of Odd(_ff„). We show in section |^ that with our definition for the FW sequence 
roughly N iterations are needed to arrive at a nonrelativistic hamiltonian that is even to order (\p\/m)^ 3 

To start the FW sequence, choose 

H = m ( 1 a +Po + O +S n ) , (9) 

where mo := m and the terms Vo, Oq, and £o are defined in Eq. This decomposition of Ho into four parts has 
the following useful properties: (i) Vq and Oq are odd; (ii) £q is even; (iii) Oq and £q are first order in parameters 
for Lorcntz violation, so products of these quantities can be neglected; and (iv) V§ is proportional to the 4x4 
Dirac identity matrix, with proportionality coefficient t\ = \p\ /m 2 . We choose for the initial FW transformation the 
hermitian matrix Sq defined by 

i5 :=^7 [Odd(ffo)] = |7°(n + Oo) ■ (10) 

This choice ensures that the odd part of exp[ad(zS'o)]i?o is smaller than the odd part of Hq. 
Our FW sequence is then defined iteratively by 

H n+1 := e iS "H n e-* S " = V 1 [»5 B> [iS n , ■ ■ ■ [iS n , H ] ■■■}}= exp[a,d(iS n )]H n (11) 

t^ k] " • ' 

k commutations with iS n 



and 



Note that 



iS n+1 :=- 7° x Odd(ff„ +1 ) . (12) 

2m n+1 



Hn+i = I J] ex P [ad(iS fe )] L H , (13) 
U=o J 

where the product represents map composition. 

In the next section, we find that each H n+ \ can be written in the form 

H n+1 = m„+i(7° + V n+ i + O n+1 + £ n+1 ) , (14) 

where the decomposition has the following useful properties: (i) V n +x and C n +i are odd; (ii) £ n +\ is even; (iii) O n +\ 
and £ n +\ are first order in parameters for Lorentz violation; and (iv) V^+i is proportional to the identity matrix, with 



proportionality coefficient determined by t\. The existence of a decomposition of the form (14) for arbitrary n, 



as well as the case (Bf) above, is a key feature making it feasible to calculate the quantum particle hamiltonian. 
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IV. CALCULATION OF THE FW SEQUENCE 



To calculate the FW sequence denned in the previous section, the explicit form is needed of the operator exp[ad(iS f „)] 
connecting H n to H n+ i according to Eq. (|i~l|). Although &d(iS n )H n can be obtained directly using the properties of 
the Dirac matrices, calculation of exp[ad(iSj)]i? n is more challenging because it is defined by an infinite series. To 
address this issue, we adopt the following approach: regard &d(iS n ) as a linear map on a suitable vector space V n 
containing both H n and H n+ i, and find a matrix expression of this map that can be exponentiated. 

The first step in implementing this approach is to define V n for each n. It is convenient to introduce V„ as the 
span of a set of basis vectors B n , defined in terms of the operators 7 , V n , O n , £ n determining H n together with the 
particular combinations of these four operators that determine ad(iS n )H n and thus also H n+ i. For each n, we define 
the ordered set 



(15) 



The eight-dimensional vector space V n is formally defined as the real span of this set, so the elements of B n by definition 
form a (linearly independent) basis. One advantage of this vector space is its relatively small dimensionality, which 
makes it susceptible to practical calculation. We can thus specify a vector V 6 V n by eight components V\ , . . . , Vg : 

V := Va Q + V 2 V„ + V 3 O n + ViVniPn, O n ] + V 5 J°[V n , £ n ] + V 6 £ n + V 7l °{V n , O n ] + V 8 V n [P n , S n ] 

"(V U ...,V S ) . (16) 

For example, H n «-> m n (l, 1, 1, 0, 0, 1, 0, 0). 

The reader is warned to avoid confusing the properties of the elements (|l^) as a basis for the vector space V n 
with their possible relationships when viewed as operators on the Hilbert space of wave functions. For example, 
the calculations below hold even if certain basis elements vanish as operators. Note also that for different n the 
corresponding vector spaces V n differ a priori. However, since both H n £ V n and H n+ \ G V„, the vector space V n is 
invariant under the action of exp[ad(iS' n )], which means V n 3 V n +i for all n. 

With the above notation, we can present the results of a direct calculation of ad(iS n )V for V G V„ performed using 
the properties of the Dirac matrices: 



zd(iS n )V «-> (t 2 n V 2 , -V u -Vi, -v 7 , \v 6 +t 2 n v & , 0, \v 2 + \v z +t 2 n v i , -V 5 ) 

In this expression, t 2 n is determined iteratively from t 2 l _ 1 through the relation 

/ cos t r . 



sin t n 



cos t n + t n sin t n 



(17) 



(18) 



Here and in what follows, we define functions of t n through their power-series expressions. All relevant functions of t n 
implicitly involve only powers of t 2 n (and hence powers of t$ = \p\ 2 /m 2 ), so it suffices to define t^. Note that t n ^ 



to leading order in i„, so t n ~ . This means that t n rapidly approaches zero if t 1, which ultimately is the 
reason for the rapid convergence of our FW sequence. 

With respect to the basis i?„, the matrix map of &d(iS„) can be extracted from Eq. (^) and is given by 



ad(iS'„) 



( 


f2 

' 1! 
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-1 























-1 









































-1 
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2 





t 2 





























1 
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1 
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t 2 














V o 











-1 








0/ 



(19) 



The exponential of this matrix can be found in closed form, but its detailed expression is unimportant. It can be used 
to calculate exp[&d(iS n )]H n , which allows us to express H n+ \ in terms of H n according to Eq. ( |l4j ) with 
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m„ + i = (cost n + t n smt n ) m n , m n+ \P n+ \ = y:ost n - — j m„:P n , 

Wn+iCri+i = I cosi„ — J m n O n + ^— t n smt n - cost„ J m n 7^ l {7 : ' n ,0 n } + -^—m n j v [V ni t n \ , 

m n+1 £ n+ i = m n £ n + (| cosi„) m„7°{7 ;, „, C„} + ( ^ j m„"P„ [Pn,£n] • (20) 

A measure of the convergence of the FW sequence can be introduced using t n ~ .In terms of a suitable matrix 
norm, ||Odd(H n+1 )|| ~ t 2 n ||Odd(# n )||+t n ||£ || ~ tf^ ||Odd(# n ) || +t ( f ] \\£ \\. Thus, as n grows ||Odd(tf„)|| rapidly 
approaches zero as (|p|/m) 3 ". Even a relatively small value of n can therefore produce a good approximation to the 
quantum particle hamiltonian. 



V. NONRELATIVISTIC QUANTUM HAMILTONIAN 

The quantum particle hamiltonian h ia \ and its nonrelativistic quantum limit h are generated in the limit of the FW 
sequence studied in the last section. Next, we demonstrate how to obtain these using simple matrix multiplication, 
and we explicitly present h rc \ and h to order t$. 

The calculation at the fcth-iteration level in the FW sequence requires obtaining the composite map 
Iln=o exp[ad(iS n )]- For each n in the FW sequence, the matrix ad(iS n ) and the action of exp[ad(iS n )] are given with 
respect to the basis B n . Since in general the vector space V„ varies with n, immediate calculation of JJ n=a exp[ad(iS' n )] 
by matrix multiplication is inappropriate. Instead, we first obtain the components of each matrix with respect to the 
special basis Bo- Ordinary matrix multiplication can then be used to derive n«=o exp[ad(i5„)]. 

The matrix for each map exp[ad(iS'„)] can be expressed in terms of t n . Explicitly, the nonzero entries for exp[ad(iS' n )] 
with respect to the basis £> are: 



exp[ad(iS„)] 



/ c„ 


+ t n s n 























\ 





























































(«„ 


— t n s n — c n )/2ig 


-t 2 s 































t n s n /2t 





t()t n S n 





















1 




















t()t n C n 








c n 







V 
















(c„ - l)/2tg 







/ 



(21) 



1 - \t 2 n and s n := sin i n /i„ 



where we have defined c„ := cost 

Since t n — > as n — > oo, it follows that exp[ad(z<S, 
this limit. The product n„ =0 exp[ad(iS„)] therefore converges as k 
quantum particle hamiltonian indeed exists. It can be shown that 



1 - if 2 
1 6 «■ 



becomes a diagonal matrix with entries (1, 0, 0, 0, 0, 1, 1, 1) in 
oo, so the limiting FW sequence giving the 



Q exp[ad(i<Sj) 



/7 











o o £ m 

\ 











-l+\t 2 ] 









1 



1 

27(7+1) 



[1 - ¥ 2 ] 











[1 - \i 



2"ol 



(22) 



where 7 := y/l + t^ is the usual relativistic gamma factor. Since the limiting FW hamiltonian is obtained by applying 
this matrix to Ho <-> (1, 1, 1, 0, 0, 1, 0, 0), the entries in brackets are irrelevant and so we have evaluated them only to 
order tjj. Thus, we find the limiting FW hamiltonian to be 



H = 7TO07 + m £o 



^i {roM 
27 



mo 



27(7+1) 



Vo[Vo,£o] 



(23) 
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an expression that is accurate to all orders in ty. Substitution from Eq. (Q) yields the explicit form 



+ i - (md j0 + ^e kl j H kl ) + 



m 



7 



+ me lm j{\gimk - Vkmgioo) 



(t — 1)771 

m(d i + d lQ ) - ^ / {mdio + \e mn iH mn ) 



P 



1 



Vjk 
P 



— mei) m( 


iffzoo) 


pk 

m 


/ k 
p p 




[(7-1)7 


7m 2 


+ 











7m 2 



-me nq ig nqk 



pkplpT: 

7m 3 



>757 



+ | [— 6j - \rne kl jg km ] - + [e l kj H m +m(dj k - r)jkd o)] ^ 

■me m i J {g m0k + g m ko) H ?7ji(^ + 2 m£m "tffmno) 

(7 - l)m 2 



p .li.pl 
7?7l 2 



P 



-m(d k i - Tjkido 



pkplptn 

1 S 
7m 13 



■757°7 J 



(24) 



This equation gives the FW form of the relativistic quantum hamiltonian for a four-component fermion. 

Certain limiting forms of the above expression are directly relevant to experiment. For applications involving 
relativistic two-component particles, such as the analysis of muon storage-ring experiments Jl4| , it suffices to retain 
only the upper left block h YC \ of H: 

U , I mC °° \ , f ! , \ ] & I \P 3 P k 

hrel = 7 m + ( fl o rne ) + dj - jm.(c j + c j0 ) - mej m{c jk - m k c 00 ) 7 

7 L J J ' Jj ryfyi J J "fm^ 



1 



-bj + mdja + \e kl .jH k i - —me Kl ^gw 

T] jk b Q + m(d jk - rijkdoo) + e l kj H m - jme lm j{\gi mk - rj km gioo) 
(7 — l)m 2 



1 



kl 



P 

7m 



-(b k + md kQ + \e mn k H mn + \me mn kg m no)Vjl 



-m(dok + d k o)r)ji + me m i ] (g m0k + g m m) 
(7 - l)m 2 r 



pkpi 
7m 2 



m(d k i - r] k id o) - \m£ nq ig nqk 



>lj, 



pkplptn 

1 3 
7m' 3 



(25) 



This is the quantum particle hamiltonian associated with the original Lorentz-violating theory. 

For many low-energy applications, including analyses of high- precision atomic experiments [[s|-[s|fi"5"|| , only nonrela- 
tivistic and subleading relativistic terms in the quantum particle hamiltonian are needed. To third order in |p|/m, 
the nonrelativistic quantum hamiltonian h for the two-component fermion is 



h = m + it— 

2m 



+ (oq - mc 00 - me ) + (-bj + mdj - \me jk ig km + \e jk iH k i) a J + [-dj + m(c 0j + c j0 ) + mej] 

+ [boSjk ~ m(d k j + d oa S jk ) - me Um {\g m ij + g m oo^ji) - £jkiH la ] ^a k + [m(-c jk - |coo<5jfc)] 
+ { [m(doj + djo) — \ (bj + mdj H m n)\ &kl 

\ me lmngmno) &jk — m£ji m (g mQk + g m ko)} 

,1, r X X PjPkPl 



m 
PjPk 



PjPk^i 



+\ [(— boSjm + mdmj + £jmnH n0 ) Ski + (-mdjk - \me knp g np j) 5w, 



PjPkPia" 
m 3 



(26) 



6 



Note that the form of Eq. (23) includes all even elements of the basis set Bq. This means that all possible combinations 
of the parameters for Lorentz violation are already contained in Eq. (^6|) . Higher-order corrections to the nonrelativistic 
hamiltonian involve only products of these combinations with powers of \p\ 2 /m 2 . One interesting implication of this 
result is that nonrelativistic experiments with single free fermions (or fermions in weak external fields) can at most be 
sensitive to the particular linear combinations of parameters for Lorentz violation appearing in Eq. ( p6| ) . Disentangling 
individual parameters requires a different class of experiment. 

As a final remark, note that our methods can also be used to obtain the nonrelativistic quantum hamiltonian h for the 
antifermion. The result for h can be expressed in the same form as Eq. (26), with the substitutions a^—t = — a M , 

0|U ► bp = +6/1; C^y > Cfiy — ~VC^ V , (1 ► d^y = G^,„ , C ^ > — ffj, ► = f fj, , QX/^V > 9 XfiV = "FffA/iJV) 

H^—y Ufa, = —H^u. This result is useful for experiments testing Lorentz symmetry with antimatter. 
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